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Abstract 

Let x%, . . . , x n be real numbers, P(x) = p n {x — xi) ■■■{x — 
x n ), and Q(x) be a polynomial of degree less than or equal 
to n. Denote by A(Q) the matrix of generalized divided 
differences of Q(x) with nodes x\,...,x n and by B(P,Q) 
the Bezout matrix (Bezoutiant) of P and Q. A relationship 
between the corresponding principal minors, counted from 
the right-hand lower corner, of the matrices B(P,Q) and 
A(Q) is established. It implies that if the principal minors 
of the matrix of divided differences of a function g(x) are 
positive or have alternating signs then the roots of the New- 
ton's interpolation polynomial of g are real and separated 
by the nodes of interpolation. 
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1 Introduction. 

In this paper a relationship between two well known matrices is es- 
tablished. The first one is a Bezout matrix B playing an important 
role in the theory of separation of polynomial roots. The second 
one is Newton's matrix of divided differences A or, in the case of 



multiple nodes, Hermite's matrix of generalized divided differences, 
playing an important role in numerical analysis and approximation 
theory. In this paper we show that the corresponding principal 
minors of B and A counted from the right-hand lower corner are 
related by a simple formula (are equal when p n — 1). An alternative 
proof of this result can be obtained from the results of [10]. As a 
simple application of the relationship between B and A, a theorem 
about locations of the roots of interpolation polynomials in terms 
of the principal minors of A is established. Many applications of 
Bezout matrix can be found in [1]-[11]. The results of this paper 
were announced without proofs in [12]. 



2 Main Results. 



With the polynomials P(x) = J2 PjX^ and Q(x) = J2 QjX-* let us 

3=0 3=0 

associate the bilinear form 



P(x)Q(y) - P(y)Q(x) 



(1) 



which Sylvester [1] named " Bezoutiant" . If the degree of Q is less 

m 

than the degree of P, that is Q(x) = 1j x \ m < n, then one adds 

3=0 

zero coefficients q m +i, . . . ,q n to Q. In what follows we assume that 
m <n and denote B(P,Q) = \ \bij\ \i j=i n . 



It has been shown (see [7]) that 
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The main properties of the Bezoutiant are (see [5, 7, 8, 9]) 



• The defect of the Bezoutiant equals the degree of the greatest 
common divisor of the polynomials P and Q. 

• The rank of the Bezoutiant matrix equals the degree of the 
last principal minor of the matrix B = ||&j,j||i,j=i,..., n which 
does not vanish if, in constructing the consecutive major mi- 
nors, one starts from the lower right-hand corner. 

• If the Bezoutiant matrix is positive definite then both poly- 
nomials P(x) and Q(x) have real, distinct roots. Moreover, 
the roots of P(x) and Q(x) interlace. 

• If all consecutive principal minors starting from the lower 
right-hand corner are positive or have alternating signs, then 
the roots of P(x) and Q(x) are real, distinct, and interlace. 

Since principal minors of Bezoutiants play so important a role, 
it seems interesting to find explicit formulas for them. If the roots 
Xi, x 2 , ■ ■ ■ , x n of P(x) are simple, such formulas were established in 



the lower right corner of the Bezoutiant B(P,Q) of polynomials 
P(x) = p n (x — x\) . . . (x — x n ) and Q(x). Then, 



[13]. 



Theorem 1 Let \ba 



ij\i,j=k+l 



be the principal minors counted from 



'i,j\i,j=k+l ~ P>, 



n 



2(n-k) 




Q(x h ) ■ ■ -Q(x in _ k ) 
P'{ Xil )---P'{x in _ k ] 



) 



X 



n 

{hi 32) c (ii, 

jl < 32 




(3) 




Remark 1 If k = n — 1 then formula (3) becomes 



Kn = Pi E 



n 



Q{xj) 

P'(XiY 



(4) 
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Remark 2 Since |^j|n = fc+i are continuous functions ofx\,..., x m 
in case of multiple roots one has to find the corresponding limit 
which is technically difficult and leads to complicated expressions. 

In order to consider the case of which are not neces- 

sarily different, let us introduce the following generalized divided 
differences. 

Definition 1 (see [14])-' 

g[x t ] := g(x t ), i = l,...,n, 



' X ik\ 



g[x 



«2 '•••»***! 



if %i\ 7^ •''ijt 



(5) 



- cia;5'[ a '' ^2 ) • • • ) •Eifc_i]x=Xj 1 j if %h 



X 



'Ik ■ 



Remark 3 This definition of generalized divided differences is equiv- 
alent to the definition given in [14] if Xi < x<i < . . . < x n . 

Consider the following triangular matrix of the generalized di- 
vided differences: A(g) = \ | Ay| \ij=i,..., n -> where 



0, if i + j < n + I 
g[x n -i+i,...,Xj], if i+j>n+ 1, 



(6) 



that is 



A($) 





A n _i 



'n- l.n 



. . 

VAi Ai, 2 . . . Ai, n _i Ai, n J 



(7) 



Remark 4 yls it is well known, Newton-Hermite's interpolation 
polynomial for n nodes {xi, . . . , x n } is A 1 + A lj2 (a: — x±) + . . . + 
Ai >n (a; - xi) . . . (x - x n _i). 

Denote by |Ajj|" 3 - =fc+1 the principal minors of the matrix A 
counted from the lower right corner. The following theorem estab- 
lishes a relationship between principal minors of the Bezoutiant and 
Newton's matrix. 
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Theorem 2 Let \bi,j\'ij=k+i an d j=k+i ^ e ^ e principal mi- 

nors of the matrices B(P,Q) and A(Q) counted from the lower 
right corner. Then 

\hJt J=k+1 = p n n - k \A hJ \l J=k+1 , k = 0, 1, . . . , n - 1. (8) 

The relationship between B(P, Q) and A(Q) established in this 
theorem is surprising taking into account that these matrices are 
of very different type, the Bezoutiant is a symmetric matrix and 
Newton's matrix is a triangular matrix. Two simple examples below 
show these matrices for some polynomials of degree three. 

Example 1 Let us consider polynomials P(x) = x 3 — Ax 2 — x + A = 
(x + l)(x — l)(x — 4) and Q(x) = x 3 — 6x 2 + llx — 6 = (x — l)(x — 

2) (x-3). Then 

( -38 48 -10 \ / 6 \ 

B(P,Q)=\ 48 -60 12 , A(Q)= 2 . 

V —10 12 -2 / V-24 12 -2 / 

Since p 3 — 1 the corresponding principal minors of these two ma- 
trices counted from the lower right-hand corner are equal, they are 
-2,-24,0. 

Example 2 Consider polynomials P(x) = x 3 — 12x 2 + 44z — 48 = 
(x-2)(x-A)(x-Q) andQ(x) = x 3 - 9x 2 + 23x - 15 = (x-l)(x- 

3) (x-5). Then 

f AAA -252 33 \ /0 15 \ 

B(P,Q) = -252 153 -21 , A(Q) = -3 9 . 

V 33 -21 3 J \3 -3 3 / 

Principal minors counted from the lower right-hand corner are 3, 18, 135. 

Theorem 2 and the properties of the Bezoutiant described above 
imply the following theorem. 

Theorem 3 If all consecutive principal minors of the matrix of 
divided differences (see (7)) of some function g(x) starting from the 
lower right-hand corner are positive or have alternating signs, then 
the roots of Newton's interpolation polynomial are real, distinct, and 
interlace with the nodes of interpolation. 
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3 Proofs. 



As it is shown in [7], 
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Since Q is a polynomial of degree m, \bij\fj =k+1 = for k = m + 
1, . . . , n — 1. Thus, let us assume that k < m and k < n — 1. Then, 



(n-k)ln-k-l) 



(9) 



where 
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This determinant can be represented as, 



d 



Pn ■ 
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■ ■ Po 
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(10) 

where Ik is k x k unit matrix (obviously, there are no rows below 
the second dashed line if k — 0). 

First let us assume that the roots of the polynomials P and Q 
are simple and distinct. 

Denote by Vj(x±, . . . , x n , yi, . . . , y n -k) the following matrix: 

/x{ ... x{ y{ ... y 3 m _ k \ 

V 1 ... 1 1 ... 1 / 

Then V n+m -k-i(xi, . . . , x n , y±, . . . , y n -k) is the Vandermont matrix 
and 

det(V r n +m-fe-i(^i, . . . ,x n ,yi, . . . , y m -k)) 

l<ii<i2<m—k 



i<h<i2<n l<h<n 

1 < %2 < m — k 



' -^\n(m—k) 

~ k f(Vl) ■ ■ ■ fiVm-k) II (Zii-Zia) II (Vii-Vh)- 

l<ii<i2<n l<ji<j2<m—k 

(11) 



rn 
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Multiplying determinants(lO) and (11) one gets: 
d ■ Vand = 



Mi 
M 2 
V fc _i(xi, ...,x n ) Vfc_i(yi, . . . ,y m ), 



where 



(yr^PiVi) ••■ vZ--t 1 P(Vm-k)\ 



m—k—l 



M 1 = 



Mo 



V P(vi) 

xr k ~ l Q{xi) 
Q(xi) 



P{y m -k) i 
x n n - k - l Q(x n ) \ 



Therefore, 



d . Vand = (-l) n ^P( yi ) . . . P(y m - k )D (y n - y j2 ), 

l<ji<j'2<m-k 

(12) 

where 



D 



x^QM ... xl~ k - l Q(x n ) 



Q(xi) 



x 



fc-1 



Q(x n ) 

„fc-l 



, if k > 1, (13) 



and 



Q(x n ) 



if jfc = o. 



(14) 



QOn) ... 

From (9), (11), and (12) one gets 

^_^(n-k)(n-k-l)/2pn-k 



■k+1 



F[ {Xii Xi 2 ^j 
l<il<i2<n 



D. 



(15) 



Since the case k = is trivial, let us assume that k > 1. Sub- 
tracting from all rows of the matrix D, except of the (n — fc)th row 
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and of the last row, the next row, multiplied by x±, and from the 
(n — k)th row the last row multiplied by Q{x\) one gets 

xr^^-x^Q^) ... x n ~ k -\x n - x x )Q{x n ) 



D 



(x 2 - x l )Q(x 2 ) 

Q(x 2 ) - Q( Xl ) 

xt 2 (x 2 - xi) 

(x 2 - xi) 

1 1 



{x n 3'l)Q(2'n) 
Q{x n ) - Q(X X ) 
X n (x n X\) 

(x n - Xi) 
1 



(16) 

After pulling out the common multipliers Xj — x±, j = 2, . . . , n from 
columns, one obtains: 

2 Q(x 2 ) ... xZ~ k ~ 2 Q(x n ) 



d = (-ir 1 n & - *o 

2<j<n 



~n-h-2, 
X 2 



Q(x 2 ) 
Q[x 1 ,x 2 ] 

x 2 



Q(Xn) 

Q[x u x n 



x 



k-2 



(17) 

Now let us consider two cases: n — k > k — 1 and n — k < k—1. 
Denote Q[i] := Q(x t ), Q[ii,i 2 , . . . , i k ] ■= Q[x h ,x i2 , . . .,x ik ]. 

The first case. If n — k > k — 1, then, continuing this process, 
after k steps one gets 

d=(-i) 2 n ( x h~ x h) 
i<ji<k 

ji <32<n 
x n k - + l k - l Q[k+l] ... xl~^Q[n] 



x 



Q[k + 1] 
Q[k,k + 1] 
Q[k-l,k,k + l] 

Q[l,...,k,k + 1] 



Q[n] 
Q[k, n] 
Q[k — 1, k, n] 

Q[l,...,k,n] 



(18) 
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Subtracting from the first n — 2k — 1 rows the next row multiplied 
by Xfc+i one obtains: 



v fc(2n-fc + 3) 

D = (-1) ^ det(C) 



n ( x j2-xn), (19) 
i < ji < fc + 1 

Ji <32<n 



where the columns of the matrix C are: 

/X n-2fc-2 

/ \ 



(x n k+ f \x k+l -x k+1 )Q[k + i]\ 



Ci 





Q[fc + 1] 
Q[fc, fc + 1] 

Vq[i,...,M + i]/ 



V 



(ajfc+i - x k+1 )Q[k + i] 

+ A 

Q[k,k + i] 
Q[k-l,k,k + i] 

Q[l,...,k,k + 2] 



J 



i = 2, . . . , n — k. (The notation C will be used below to denote 
different matrices.) 

After subtracting the first column from the other columns and 
pulling out the common factors (x k +i — x k ), . . . , (x n — x k ), one gets 



fc(2n-fc + 3) 



n 

i < ji < k + 1 

ji <32<n 



x 



3 2 



x jl ) x 







^; 2 2fe " 2 Q[^ + 2] 



X 



n-2fc-2 



Q[n] 



Q[fc + 2] 

Q[fc + 1] Q[fc + l,fc + 2] 

g[fc,fc + l] Q[fc,fc + 2] 

Q[Jfc-l,fc,fc+l] Q[k- l,fc,fc + 2] 

Q[l,...,fc,fc + 1] Q[l,...,fc,fc + 2] 



Q[fc + l,ra] 
— 1, fc, n] 

Q[l, • • • , fc,n] 



(20) 

Again, let us subtract from the first n — 2k — 2 rows the next row 
multiplied by x k+ 2, then, subtract the second column from the other 
columns, and pull out from the columns common factors (x k+ 2 ~ 
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Xk+i), • • • , (x n — Xfc+i). Continuing this process, after n — 2k — 1 
steps one gets: 



, fc(2n-fe + 3) 

D = (-1) ^ det(C) 



n i x j2 x ji)i 

l<ji<n — k — 1 
ji <32<n 



where the columns of C are 
/ 



d = 





Q[k + i] 
Q[k + i-l,k + i] 

V Q[l,...,A; + i] J 

( Q[k + 1] \ 

Q[n — k — 1, k + % 



i — 1, . . . , n — 2k — 1, 



, i = n — 2k, . . . ,n — k. 



\Q[l,...,n-k- l,k + i]J 
Let us subtract from the last k columns the previous column and 

n 

factor out J] ( x j ~ x j-i). Then let us repeat this procedure 

j=n— k+1 

with the last k — 1 columns, and so on. Finally, 

fc(2n-fc + 3) „ 

^ = (-1) 2 II ( X j2-X h )x 

l<jl<j2<m 






Q[n- 




• • • Q[-^n— fej • 


■ ■ 1 x n 




Q[x k ,x k+ i] 


.. Q[x k+1 ,... 

■ ■ Q[ x ki ■ ■ ■ 


, x n — k g 

x n—k\ 


... Q[x k +i,. 
■ ■ ■ Q[x k , ■ ■ 


■ ■ i Xn 

■ ■ i Xn 

■ i x n ] 


Q[x u . . . ,x k +i] ■ 


.. Q[X!,... 


Xn— k] 


Q[X!,.. 


■ -i x n \ 


Since, n ( x j 2 


~ x n) = (-1) 


n(n-l)/2 


n {x h - 

l<jl<j2<n 


X 32 ) ) 



(21) 



obtains (8) from (6), (15) and (21). 
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The second case. If n — k < k — 1, then, transforming the 
determinant similarly to the first case, from (17) one obtains: 



tn+k+4)(n-k-l) 

D =(-!)- 2 1 det(C) 



where the columns of C are 
/ 



n 

l<ji<n — k — 1 
ji <j2<n 



(x 



32 



Xj 1 ) , 



Q[n — k — l + i] 
Q[l,n- k — l + i] 

Q[l, . . . , n — k — 1, n — k — 1 + i] 



x. 



2k-n 
n— k—l+i 



V 



•En— k—l+i 
1 



l,...,k+ 1. 



/ 



After subtracting from the first n — k rows the last row multiplied 
by the first element of the row, from the n — k + 1th to fcth rows the 
next row, multiplied by and pulling out the common factors 

•En—k+l "En—kf • • i -^n -^n— k One gets 



v (n+fc+3)(n-fe) , v 

D = (-1)- ^ — i det(C) 



n i x j2 x ji)i 

1 < j 'i < n — k 
ji <32<n 



d = 



( Q[n — k,n — k + i] \ 
Q[n — k — l,n — k,n — k + i] 



Q[l, . . . ,n — k,n — k + i] 

2k-n-l 



x 



V 



n—k+i 



%n—k+i 
1 



, i — 1, . . . , k. 



Continuing this process, after 2k — n steps one obtains: 



(fe-l)(2n-fe+4) 

D = (-l)^-^ i det(C) 



n ( x 32 x ji , 

l<ji<k-l 
ji <32<n 



12 



( Q[n-k,...,k-l,k-l + i] \ 
Q[n - k - 1, . . . , k - 1, k - 1 + i\ 



V 



Q[l,...,k-l,k-l + i] 
1 



, i — 1, . . . , n — k + 1. 



After subtracting the first column from the next columns and pulling 
out the common factors x^+i — £&,. . . , x n — Xk one gets 



v fc(2n-fc + 3) , v 

D = (-1) 2 det(C) 



l<ji<k 
ji <32<n 



( Q[n — k, . . . ,k,k + i] \ 
Q[n — k — 1, n — k, . . . , k, k + i 



V 



, i — 1, . . . , n — k. 



J 



Q[l,...,k,k + i] 

Finally, one has to subtract the first column from the next columns 
and pull out the common factors Xk+2 — • • , x n — Xk+i, then 

one has to subtract the second column from the next columns and 
pull out the common factors Xk+z — Xk+2, - ■ ■ , x n — and so on. 
Then one gets 



fc(2n-fc + 3) 

D = (-l)^— 



n ( x j2 x ji ) 

1 < ji < 32 < n 



x 



^n-k,k+l{Q) ^n-k,k+2{Q) ■ ■ ■ ^n-k,n{Q) 

(Q) 



A 1>k+1 (Q) A lifc+2 (Q) 
Then, (15) and (22) imply (8). 



A 1>n (Q) 



(22) 



Thus, Theorem 2 has proven in the case of simple and distinct 
roots of polynomials P and Q. 

To prove the theorem in the general situation, we show first that 
it remains true if polynomials have common roots. 

n r m—r 

Let P(x) = p n n (x - xj) and Q(x) = q m U (x - xj) U (x - 
j=i j=i j=i 

Uj) for some r, < r < min(m, n). Then for sufficiently small e 
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polynomials P(x) and Q e {x) = q m I\ (x — Xj — e) II (x — Uj) have 

3=1 3=1 

distinct roots and therefore (8) holds. Since both sides in (8) are 
continuous functions of e, the formula remains true when e — > 0. 

Similarly one can prove (8) in the case of multiple roots. Denote 
by r the highest multiplicity of the roots of P(x). We will use 
induction with respect to r. If r = 1 the roots of P(x) are simple. 
Assume that (8) is true for some r and prove it for r + 1. Assume 
that there is one root of multiplicity r + 1. (The case of several 
roots of multiplicity r + 1 can be proved similarly.) Let P(x) = 
(x—Xi) r+1 Pi(x), where P\{x) is a polynomial of degree n— r— 1 with 
the roots distinct from x±. Let P e (x) = (x — x\ — e)(x — xi) r Pi(x) 
for sufficiently small e. By assumption, (8) is true for P £ and Q. 
The left hand side in (8) is a continuous function of e. In the right 
hand side 

lima[zi + e,xx, . . . , x 1 , x 2 , . . . , x a ] = g[x lt ...,x 1 ,x 2 ,.. .,x a \. 

£->0 v v ' ^ v ' 

r r+1 

This observation completes the proof of Theorem 2 in the case of 
multiple roots. 



References 

[1] J. Sylvester, On a Theory of the Syzygetic relations of two 
rational integral functions, comprising an application to the 
theory of Sturm's Functions, and that of the greatest Alge- 
braic Common Measure, Philos. Trans. Roy. Soc. London 143 
(1853), 407-548. 

[2] C. Hermite, Extrait d'une lettre de Mr. Ch. Hermite de Paris 
a Mr. Borchardt de Berlin, sur le nombre des racines d'une 
equation algebrique comprises entre des limites donnees, J. 
Reine Angew. Math. 52 (1856), 39-51. 

[3] A. Cayley, Note sur la methode d 'elimination de Bezout, J. 
Reine Angew. Math. 53 (1857), 366-376. 

[4] A. Hurwitz, Ueber die Bedingungen unter welchen eine Gle- 
ichung nur Wurzeln mit negativen reellen Teilen besitzt, Math. 
Ann. 46 (1895), 273-284. 



14 



[5] M.G. Krein and M.A. Naimark, The method of symmetric and 
Hermitian forms in the theory of the separation of the roots of 
algebraic equations, Linear and multilinear algebra, 10 (1981), 
265-308 (The paper was originally published in Kharkov in 
1936). 

[6] A.S. Householder, Bezoutiants, elimination and localization, 
SIAM Review, 12 (1970), NO. 1, 73-78. 

[7] F. I. Lander, Bezoutiant and inversion of hankel and toplitz 
matrices, Matematicheskie Issledovaniya, 9, N2 (32), (1974), 
69-87. 

[8] P.A. Fahrmann and B.N. Datta, On Bezoutians, Van der 
Monde matrices, and the Lienard-Chipart stability criterion, 
Linear Algebra and its applications, 120 (1989), 23-37. 

[9] A. Olshevsky and V. Olshevsky, Kharitonovs theorem and Be- 
zoutians, Linear Algebra Appl., 399, (2005), 285297. 

[10] G.M. Diaz-Toca, L. Gonzalez- Vega, Various New Expressions 
for Subresultants and Their Applications, Applicable Alge- 
bra in Engineering, communication and Computing 15 (2004), 
233-266. 

[11] C. D'Andrea, H. Hong, T. Krick and A. Szanto, An elementary 
proof of Sylvester's double sums for subresultants, Journal of 
Symbolic Computation, 42 (2007), 290-207. 

[12] R.G. Airapetyan, The relationship between Bezoutian matrix 
and Newton's matrix of divided differences, Proceedings of the 
7th International ISAAC Congress, Imperial College, London, 
UK, 2010, 573-578. 

[13] R.G. Airapetyan, On the reduction of the Cauchy problem for 
a hyperbolic equation to symmetric systems, Soviet Journal of 
Contemporary Mathematical Analysis, 21 (1986), N 1. 

[14] J. Stoer and R. Bulirsch, Introduction to Numerical Analysis, 
3rd ed., Texts in Applied Mathematics, Springer, New York, 
2002. 



15 



